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The electromagnetic radiation by quark-gluon plasma in strong magnetic field is calculated. 
The contributing processes are synchrotron radiation and one~photon annihilation. It is 
shown that in relativistic heavy-ion collisions at RHIC and LHC synchrotron radiation 
dominates over the annihilation. Moreover, it constitutes a significant part of all photons 
produced by the plasma at low transverse momenta; its magnitude depends on the plasma 
temperature and the magnetic field strength. Electromagnetic radiation in magnetic field 
is probably the missing piece that resolves a discrepancy between the theoretical models 
and the experimental data. It is argued that electromagnetic radiation increases with the 
magnetic field strength and plasma temperature. 



I. INTRODUCTION 



It has been known for a long time that strong magnetic fields can be generated in heavy-ion 
collisions. Recent calculations [IHS] confirm this wisdom and predict, using essentially classical 
electrodynamics, that the magnetic field shortly after the collision reaches 10^® — 10^^ G, which 
by far exceeds the critical value Be = 4.41 • 10^^ G for electrons, known also as the Schwinger 
field. However, it was only recently appreciated that such fields may have great phenomenological 
significance. There are two observations leading to this conclusion: (i) phenomenological models 
suggest that the Quark Gluon Plasma (QGP) is formed after a very short time after the heavy-ion 
collision - on the order of a few tenth of a fm/c; (ii) the relaxation time of magnetic field in the 
presence of QGP is proportional to the plasma electric conductivity due to the induction of electric 
Foucault currents by the time-dependent magnetic field [6j. The relaxation time of magnetic field 
is estimated to be about 1-2 few fm/c Therefore, magnetic field has a profound influence on 
all aspects of the physics of relativistic heavy-ion collisons. In particular, it was argued in [6] that 
magnetic field induces energy loss by fast quarks and charged leptons via the synchrotron radiation 
and polarization of the fermion spectra. It contributes to enhancement of dilepton production at 
low invariant masses [7J and enhances the azimuthal anisotropy of the quark-gluon plasma (QGP) 
[HI [9]. It causes dissociation of the bound states, particularly charmonia, via ionization [10^ 111]. 
Additionally, magnetic field drives the Chiral Magnetic Effect (CME) [U [T^HT^| . which is the 
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generation of an electric field parallel to the magnetic one via the axial anomaly in the hot nuclear 
matter. Additionally, the effect of the magnetic field on the QCD phase diagram was studied using 
model calculations [TBH3H] and lattice simulations [55ti57] . It is also argued in Appendix, that at 
RHIC, at early times after heavy ion collision, about 3% of energy density of plasma resides in the 
magnetic field, while at LHC, this fraction reaches as much as 40%. 

This paper addresses the problem of photon radiation by quarks and antiquaries of QGP moving 
in external magnetic field. This radiation originates from two sources: (i) synchrotron radiation 
and (ii) quark and antiquark annihilation. QGP is transparent to the emitted electromagnetic 
radiation because its absorption coefficient is suppressed by a^. Thus, QGP is shinning in magnetic 
field. The main goal of this paper is to calculate the spectrum and angular distribution of this 
radiation. In strong magnetic field it is essential to account for quantization of fermion spectra. 
Indeed, spacing between the Landau levels is of the order \J eB, while their thermal width is of the 
order T. Spectrum quantization is negligible only if eB <^ which is barely the case at RHIC 
and certainly not the case at LHC (at least during the first few fm's of the evolution). Fermion 
spectrum quantization is important not only for hard and electromagnetic probes but also for the 
bulk properties of QGP. 

The presentation is structured as follows. In Sec. [H] the spectrum and angular distribution of 
synchrotron radiation by QGP is calculated in the ideal gas approximation and compared to the 
experimental data. The results exhibited in Fig. 3j4 indicate that photon radiation in magnetic 



field gives a significant contribution to the total photon yield. This contribution seems to be 
large enough to account for the discrepancy between the model calculations assuming no magnetic 
field and the experimental data at RHIC |38j.* Moreover, the electromagnetic radiation rapidly 



increases with B and T suggesting that it plays even more important role at LHC. In Sec. HI 
the photon spectrum emitted in pair annihilation is calculated and it is shown that it is small 
as compared to synchrotron contribution. A possible way to ascertain existence of synchrotron 
radiation is discussed in Sec. IIVI 

II. SYNCHROTRON RADIATION 

Motion of charged fermions in external magnetic field, which we will approximately treat as 
spatially homogeneous, is quasi-classical in the field direction and quantized in the reaction plane, 



* Another possible explanation has been recently suggested in [39] . 
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which is perpendicular to the magnetic field and span by the impact parameter and the heavy 
ion collision axis. In high energy physics one usually distinguishes the transverse plane, which is 
perpendicular to the collision axis and span by the magnetic field and the impact parameter. The 
notation is adopted in which three-vectors are discriminated by the bold face and their component 
along the field direction by the plain face. Momentum projections onto the transverse plane are 
denoted by subscript _L. 

In the configuration space, charged fermions move along spiral trajectories with the symmetry 
axis aligned with the field direction. Synchrotron radiation is a process of photon 7 radiation by a 
fermion / with electric charge ej = ZfC in external magnetic field B: 

f{ef,j,p)^ f{ef,k,q)+j{k), (1) 

where k is the photon momentum, p, q are the momentum components along the magnetic field 
direction and indicies j, A; = 0, 1, 2, . . . label the discrete Landau levels in the reaction plane. The 
Landau levels are given by 



Ej = yjm? + 2jefB , ek = \Jnn? + q^ + 2ke fB , (2) 

In the constant magnetic field only momentum component along the field direction is conserved. 
Thus, the conservation laws for synchrotron radiation read 

Ej = uj + Ek , p = q + u) cos 9 , (3) 

where cj is the photon energy and 9 is the photon emission angle with respect to the magnetic field. 
Intensity of the synchrotron radiation was derived in [lO] . In |41H44j it was thoroughly investigated 
as a possible mechanism for 7-ray bursts. In particular, synchrotron radiation in electromagnetic 
plasmas was calculated. Spectral intensity of angular distribution of synchrotron radiation by a 
fermion in the j'th Landau state is given by 

^ = E E {I-^^I^ + l-^ii 1^} ^i- - + (4) 

where Fj^ = (1 + Sjo){l + 5ko) accounts for the double degeneration of all Landau levels except 
the ground one. The squares of matrix elements Ai, which appear in Q, corresponding to photon 
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polarization perpendicular and parallel to the magnetic field are given by, respectively, 

4^£j£k\M±\'^ ={ejek-pq - m^)[llk-i + ^|-i,fc] + 2^2jefB^2kefB[Ij^k^iIj_i^k] ■ (5) 
4ejefe|X|,|2 =cos^ e{{ej£k - pq - rr?)[ll^_^+ - 2^2jefB^2kefB[Ij^k-iIj^i,k]} 

- 2cos9sm9{p^/2kefB[Ij-l^kIj-l,k-l + Ij,k^ilj,k] 

+ Qy^2jefB[Ij^kIj-i,k + Ij~i,k^ilj,k-i]} 

+ sin2 9{{e,ek + pq - m^)[l]_,^k-i + Ij,k] + 2^2jefBy/2kefB{Ij^i^k^iIj,k)} , (6) 

where for j > k, 

^ hki^) = (-l)^-'=y|"e-f x^Lf^x). (7) 

and Ij^k{x) = Ikj{x) when k > j. {Ij-i are identically zero). The functions L-j^^{x) are the 
generalized Laguerre polynomials. Their argument is 



2 

sin^ 9 . (8) 



a; . 2 



2efB 

Angular distribution of radiation is obtained by integrating over the photon energies and re- 
membering that Ek also depends on lo by virtue of Q and ([s]): 

^^=y:^E — . r,,{i>fjHiA^,n. (9) 

/ fe=0 ■' 

where photon energy cj is fixed to be 

u* = ^^l{ej-pcos9)-\{£j-pcos9f-2efB{j-k)sm'^9Y^H . (10) 
sin 9 J 

In the context of heavy-ion collisions the relevant observable is the differential photon spectrum. 
For ideal plasma in equilibrium each quark flavor gives the following contribution to the photon 
spectrum: 

f °° j=0 k=0 

where 2Nc accounts for quarks and antiquarks each of Nc possible colors, and (2 — 5jfi) sums over 
the initial quark spin. Index / indicates different quark flavors. V stands for the plasma volume. 
The statistical factor /(e) is 
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The 5-function appearing in Q can be re-written using ^ and Q as 

6{UJ - Ej + Ek) = , (13) 



where 



p± =■[ cos 9{m'j — m| + cj^ sin^ i 



(m-j + mfc)2 — sin^ 6] [{rrij — rrikY — uJ^ sin^ ^] | / (2w sin^ B) . (14) 
The following convenient notation was introduced: 

= 'm'^ + 2jefB , ml = m'^ + 2ke j B . (15) 



The physical meaning of ( 14 ) is that synchrotron radiation of a photon with energy uj at angle 9 by 
a fermion undergoing transition from j'th to fc'th Landau level is possible only if the initial quark 
momentum along the field direction equals p^. 

Another consequence of the conservation laws ([s]) is that for a given j and k the photon energy 



cannot exceed a certain maximal value that will be denoted by ujgjk- Indeed, inspection of (14) 
reveals that this equation has a real solution only in two cases 

(i) rrij — nik > uj sin 6 , or (ii) rrij + m/j < uj sin 6 . (16) 

The first case is relevant for the synchrotron radiation while the second one for the one-photon pair 
annihilation as discussed in the next section. Accordingly, allowed photon energies in the j — ?■ 
transition satisfy 



m 



ruk ^Jm? + 2jefB - y^m^ + 2kefB 



UJ < UJs,jk = \ . = ■ n — ■ (17) 

No synchrotron radiation is possible for uj > uj^ j/^. In particular, when j = k, UJgjk = 0, i.e. no 
photon is emitted, which is also evident in ( |10[ ). The reason is clearly seen in the frame where 
p = 0: since ej > Sk, constraints ^ and ^ hold only if a; = 0. 

Substitution of Q into ( |11[ ) yields the spectral distribution of the synchrotron radiation rate 
per unit volume 



2^ — ^^^efB2^2^uj{l + 6ko)^{uJs,ij-uj) J dp2_^ 



/ j=0 fc = ± I Efc I 

x{|A^^|2 + |A^|||2}/(e^.)[l _/(£,)], (18) 

where is the step-function. 
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The natural variables to study the synchrotron radiation are the photon energy oj and its emis- 
sion angle 9 with respect to the magnetic field. However, in high energy physics particle spectra are 
traditionally presented in terms of rapidity y (which for photons is equivalent to pseudo-rapidity) 
and transverse momentum k±. k± is a projection of three-momentum k onto the transverse plane. 
These variables are not convenient to study electromagnetic processes in external magnetic field. 
In particular, they conceal the azimuthal symmetry with respect to the magnetic field direction. 
To change variables, let z be the collision axis and y be the direction of the magnetic field. In 
spherical coordinates photon momentum is given by k = Ljj{smacos(j)x + smasiiKpy + cosaz), 
where a and (f) are the polar and azimuthal angles with respect to z-axis. The plane xz is the 
reaction plane. By definition, k ■ y = cos 9 implying that cos 9 = sin a sin (p. Thus, 



B + k?. = ^. y 



y 



In tan 



a 



sm ( 



The second of these equations is the definition of (pseudo)-rapidity. Inverting ( 19 ) yields 

sin 6 



k± cosh y , cos 9 



coshy 

Because dy = dk^/uj the photon multiplicity in a unit volume per unit time reads 



^^synch 



^^synch 



UJ 



dVdtd'^k±dy dVdtd^k dVdtLodudn 



(19) 



(20) 



(21) 
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FIG. 1: Spectrum of synchrotron radiation by u quarks at eB — rn^, y — 0, (p — 7r/3: (a) contribution of 
10 lowest Landau levels j < 10; several cutoff frequencies are indicated; (b) summed over all Landau levels. 
m„ = 3 MeV, T ^ 200 MeV. 

Fig. [l] displays the spectrum of synchrotron radiation by u quarks as a function of k± at fixed 



(j). At midrapidity y = (20) implies that k± = to. Contribution of d and s quarks is qualitatively 
similar. At eB ^ mP, quark masses do not affect the spectrum much. The main difference 
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stems from the difference in electric charge. In panel (a) only the contributions of the first ten 
Landau levels are displayed. The cutoff frequencies cosjk can be clearly seen and some of them are 
indicated on the plot for convenience. The azimuthal distribution is shown in Fig.[2j Note, that at 
midrapidity (j) = tt /2 — 0. Therefore, the figure indicates that photon production in the direction 
of magnetic field (at = 7r/2) is suppressed. More photons are produced in the direction of the 
reaction plane (p = ^. This results in the ellipticity of the photon spectrum that translates into the 
positive "elliptic flow" coefficient V2- It should be noted, that the classical synchrotron radiation 
has a similar angular distribution. 




12 6 4 3 ~f2 2 ~4 





FIG. 2: Azimuthal distribution of synchrotron radiation by u-quarks at /c_l — 0.2 GeV, eB — toJ, y = 0. 
m„ = 3 MeV. 



In order to compare the photon spectrum produced by synchrotron radiation to the photon 
spectrum measured in heavy-ion collisions, the d and s quarks contributions were summed up. 
Furthermore, the experimental data from [38] was divided by Vt, where t is the magnetic field 
relaxation time. The volume of the plasma can be estimated as F = 7ri?^t with i? ~ 5 fm being 
the nuclear radius. Therefore, 

dA^elp ^ ^iVelp 1 ^ rfA^clp /GeVy / Ifm 
dVdtd'^k±dy d'^k±dy ttRH^ d^k±dy \U.9 J \ t 

The results are plotted in Fig.jsj In panel (a) it is seen that synchrotron radiation gives a significant 
contribution to the photon production in heavy-ion collisions at RHIC energy. This contribution is 
larger at small transverse momenta. This may explain enhancement of photon production observed 
in |38j . Panel (b) indicates the increase of the photon spectrum produced by the synchrotron 
radiation mechanism at the LHC energy. This increase is due to enhancement of the magnetic field 
strength, but mostly because of increase of plasma temperature. This qualitative features can be 
better understood by considering the limiting cases of low and high photon energies. 
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(a) 



FIG. 3: Azimuthal average of the synchrotron radiation spectrum of u,d,s quarks and their corresponding 
antiquarks. (a) eB — m^, y = compared to the experimental data from [38 divided by Vt = 25tt fm^ 
(dots) and Vt^9x 25tt fm^ (stars), (b) eB = ml,T ^ 200 MeV, y = (sohd hne) compared to eB ^ 15ml, 
T = 400 MeV, y = (dashed hne). m„ = 3 MeV, rrirf = 5 MeV, m, = 92 MeV. 

A. Low photon energy 



The low energy part of the photon spectrum satisfies the condition oj <^ y^CfB. The corre- 
sponding initial quark momentum component along the field p and energy ej follow from ( 14 ) and 
([2]) and are given by 

(j - k)efB{cose±l) 



P± 



oj sm 



+ 0{lo), e,^\pl\+0{u:). 



(23) 



Evidently, Sj ^ eB. In practice, the magnetic field strength satisfies VelS > T, so that Sj ^ T. 
Therefore, synchrotron radiation is dominated by fermion transitions from low Landau levels due 



to the statistical factors appearing in (11). 



For a qualitative discussion it is sufficient to consider the 1 — )• transition. In this case the 
matrix elements ^ and ([g]) read 



\M 



1,0|2 



7r~ — |-^io(^i^o -pqcos^ 6 - w?) + cosds\ii6qJ2efBIifilQfl\ 



(24) 



Assuming that the field strength is supercritical, i.e. CfB ^ m^, but keeping all powers of oj (for 



future reference) ( 14 ) reduces to 

1 



2oj sin^ e 



{2e/B(cos6'±l)+a;^sin2(9(cos6'=Fl)} . 



(25) 
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Furthermore, using the conservation laws ^ we obtain in this approximation 

1 



2 „;„2 

1 y t-i /■ / "A I / -1 n_: n — i — ■ i / . i 

2a; sin^ 6^ 

=- — ^^(2e/5-a;^sin^6l)(cos6'±l), 



2e/-B(cos 6* ± 1) - sin^ 6'(cos 6* =F 1) 

2 „• 2 



2uj sin^ 6* 

The values of the non-vanishing matrix elements Ij^k defined by ([T]) are 

Ii,o(x) = -xi/V-/2, /o,o(x) = e--/2. 

Y^2eji?/ sin 0. Then iS^ implies x 



For j = 1, k = we write using (17) co. 



s,10 



(26) 
(27) 
(28) 

(29) 



Substituting (l25|-(l29j) into gives 

1 cos 6^(1 + x) ± (1 



\M 



1,0|2 
± I 



1 



X 2 
COS I 



2(1 -x)cos6'sin2 6l 
cos 61(1 -x)±{l + x) 



cos9{l-x)±{l + x 

According to (18) the contribution of the 1 — ?• transition to the synchrotron radiation reads 



(30) 



VdtdQduj 



.^INczla efB sr^ ^ 

E^r^-r^E/(-i)[i-/(-o)]|A^i 



,0|2 



(1 -x)cos6'± (1 + x) , 

X — ^ '&{uJs,10 - W) . 



(31) 



-2x(cos6l^ 1) 

Consider radiation spectrum at ^ = 7r/2, i.e. perpendicular to the magnetic field. The spectrum 
increases with x and reaches maximum at x = 1. Since x = uj'^/{2efB), spectrum decreases with 
increase of B at fixed uj. This feature holds at low x part of the spectrum for other emission angles 
and even for transitions form higher excited states. However, at high energies, it is no longer 
possible to approximate the spectrum by the contribution of a few low Landau levels. In that 
case the typical values of quantum numbers are j,k ^ 1. For example, to achieve the numerical 
accuracy of 5%, sum over j must run up to a certain jmax- Some values of jmax are listed in Table |l} 



B. High photon energy 

The high energy tail of the photon spectrum is quasi-classical and approximately continuous. In 
this case the Laguerre polynomials can be approximated by the Airy functions or the corresponding 
modified Bessel functions. The angular distribution of the spectrum can be found in [41j : 



^^synch zjanfojm^ jefBTsiuu 



Vdtdndu ^ TT 4T3 V m3 
/ 



e-"^/^ , (32) 
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TABLE I: The upper summation limit in ( 18 ) that yields the 5% accuracy, jmax is the highest Landau level 



of the initial quark that is taken into account at this accuracy. Throughout the table y — Q. 



provided that u ^ m^JmT j e^B sin Q. Here nj is number density of flavor /, which is independent 



Uf 



2 • INcCfB 



E 



dpe 



AN, 



(33) 



Here summation over j was replaced by integration. It follows that this part of the spectrum 
increases with magnetic field strength as \fB and and with temperature as \fTe^^l'^ . Therefore, 
variation of the spectrum with T is much stronger than with B. The T dependence is shown in 
Fig.H 



0.001 
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FIG. 4: Variation of the synchrotron spectrum with plasma temperature. Lower line: T ~ 200 MeV, upper 
line: T — 250 MeV. Other parameters are the same as in Fig. [sj^a). 



III. PAIR ANNIHILATION 

The theory of one-photon pair annihilation was developed in |45t I46j . It was shown in 07] that 
in the super-critical regime eB ^ m? one-photon annihilations is much larger than the two-photon 
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annihilation. In this section the one-photon annihilation of q and q pairs in the QGP is calculated. 
For qq pair annihilation the conservation of energy and momentum is given by 

ej + £k = ^ ■, p + q = CO cos 9 . (34) 

The spectral density of the annihilation rate per unit volume reads 

dArannih azjuN^^^ r 2efB f 2efB 

f J j=0 k=0 

x5{p + q - u cos9)5{ej + £k - ^^){\Ti_\'^ + |7^P}, (35) 

where the matrix elements T can be obtained from ([sJljQ by making substitutions £k — >• — Sfc, 
q ^ —q and are given by 

4ejefc|r±|2 ={ej£k -pq + m?)[l],k~i + Ij-i,k] " 2y'2jefBy/2kefB[Ij^k-iIj-i,k] ■ (36) 
4ej£fc|7^|2 =cos20{(ejefc + + + 2v/2ie/Sv^2fce/S[/j-fc_i/j_i,fc]} 

- 2cos^sin6'{ - p^/2kejB[Lj^i^kIj-i,k-i + Ij,k-ilj,k] 

+ Q\/yefB[Ij,kIj-i,k + Ij~i,k-ilj,k-i]} 

+ sin2 9{{£jek+pq + m^)[lf_^ ,,_^ + /|J - 2./WfBy/2kefB{Ij^i^k-iIj,k)} , (37) 

with the same functions /jj as in ([T]). Integration over q removes the delta function responsible for 
the conservation of momentum along the field direction. The remaining delta function is responsible 



for energy conservation and can be written in exactly the same form as in ( 13 ) with particle energies 



and momenta now obeying the conservation laws (34). It is straightforward to see that momentum 



is still given by (14),(15). The photon spectrum produced by annihilation of quark in state j 



with antiquark in state k has a threshold Wa,ij that is given by the case (ii) in (16): 

rrij + mk ^m? + 2jefB + ^/m?~T2keJB 

OJ > UJa,ij = = ■ [Ati) 

sm9 sm9 
Thus, the spectral density of the annihilation rate per unit volume is 



azju^ r 5{p 

f j=0 k=0 ± I Eft I 



X 



{\n\' + \rii\'}f{s,)fiek). (39) 



Passing to y and p± variables in place of u and 9 is similar to (21 ). 

The results of the numerical calculations are represented in Fig.jsj Panel (a) shows the spectrum 
of photons radiated in annihilation of u and u. We conclude that contribution of the annihilation 
channel is negligible as compared to the synchrotron radiation. 
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(a) 



FIG 

at 6 



5: Photon spectrum in one-photon annihilation of u and u quarks. eB — rn^, y — 0. (a) fcj^-spectrum 
~ 7r/3, (b) aziniuthal angule distribution at k± = I GeV. 



IV. CONCLUSIONS 



Results of the calculations performed in this article indicate that photon production by QGP due 
to its interaction with external magnetic field give a considerable contribution to the total photon 
multiplicity in heavy-ion collisions. This is seen in Fig. |3] were the model calculation is compared 
with the experimental data [38]. The two processes were considered: synchrotron radiation and pair 
annihilation. In the kinematic region relevant for the current high energy heavy-ion experiments, 
contribution of the synchrotron radiation is about two orders of magnitude larger than that of pair 
annihilation. The largest contribution to the photon multiplicity arises from photon momenta of 
the order of V eB. This may provide an explanation of the photon excess observed by the PHENIX 
experiment [38] . Similar mechanism is also responsible for enhancement of low mass di-lepton 
production that proceeds via emission of virtual photon which subsequently decays into di-lepton 
pair. 

One possible way to ascertain the contribution of electromagnetic radiation in external magnetic 
field is to isolate the azimuthally symmetric component with respect to the direction of the magnetic 
field. It seems that synchrotron radiation dominates the photon spectrum at low k±. Thus, 
azimuthal symmetry can be easily checked by simply plotting the multiplicity vs to, 9 and ip, where 
uj is photon energy, is emission angle with respect to the magnetic field and ip is azimuthal angle 
around the magnetic field direction (which is perpendicular both to the collision axis and to the 
impact parameter). In Fig. [l]^a) it is also seen that in these variables it may be possible to discern 
the cutoff frequencies Ugjk that appear as resonances (in Fig. [l] y = so = u). Note that 
averaging over the azimuthal angle a around the collision axis direction distroys these features, see 
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Fig. [3} 

The greatest source of uncertainty in the present calculation stems from treating the magnetic 
field as constant. It is inevitable that it has spatial [S] and temporal variations, which will modify 
the photon spectrum. Analytical calculations of these effects present a serious challenge, but may 
be tackled in the quasi-classical approximation. Novel computational techniques, such as discussed 
in [IH] is another promising avenue for investigating the particle production in external fields. 
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Appendix A: Magnetic component of the QGP energy density 

Energy density associated with magnetic field is (in Gauss units) 

eu = {eBf/{87ra) . (Al) 

Assume that the effect of magnetic field on the QGP is weak and its energy density is much smaller 
than the total QGP energy density e^/ <^ eqcp- Then energy density of ideal QGP containing Nf 
quark flavors at temperature T is given by 

eQGP = (2(Ar2 _ 1) + Un.Nj^ . (A2) 

With Nc = 3, Nf = 2 we arrive at the following ratio 

^ = O.45M! . (A3) 
eqcp -t 

At RHIC eB ~ and T ~ 2m^, so that at early times about 3% of energy density of plasma 
resides in the magnetic field. At LHC, eB ~ 15m^ and T ~ Am.^ so that as much as 40% is stored 
in magnetic field! This signals that magnetic field plays a crucial role in QGP dynamics at LHC 
energies. 
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